Critical exponents of random XX and XY chains: 
Exact results via random walks 



in 

(N 

c 

C/5 



C3 



I 

C 

o 
o 



> 

00 

m 

o 

On 

^t— » 
c3 



c 
o 
o 



X 
S3 



Hciko Ricger 1,2 ' 3 , Robert Juhasz 4 and Ferenc Igloi 3,4 ' 5 

1 NIC c/o Forschungszentrum Jiilich, 52425 Jiilich, Germany 

2 Institut fur Theoretische Physik, Universitat zu Koln, 50937 Koln, Germany 

3 Laboratoire de Physique des Materiaux, Universite Henri Poincare, BP 239, F-54506 Vandosvre les Nancy Cedex, France 

4 Institute for Theoretical Physics, Szeged University, H-6720 Szeged, Hungary 

5 Research Institute for Solid State Physics and Optics, H-1525 Budapest, P.O.Box 49, Hungary 

(June 14, 1999) 

We study random XY and (dimerized) XX spin-1/2 quantum spin chains at their quantum phase 
transition driven by the anisotropy and dimerization, respectively. Using exact expressions for 
magnetization, correlation functions and energy gap, obtained by the free fermion technique, the 
critical and off-critical (Griffiths-McCoy) singularities are related to persistence properties of random 
walks. In this way we determine exactly the decay exponents for surface and bulk transverse and 
longitudinal correlations, correlation length exponent and dynamical exponent. 



Disordered quantum spin chains have gained much in- 
terest recently p]-P]. It seems to be established right 
now that the critical properties in these one-dimensional 
system are governed by an infinite-disorder fixed-point 
[|9| and the application of a renormalization group (RG) 
scheme a la Dasgupta and Ma pi is a powerful tool 
to determine critical properties and static correlations of 
these new universality classes, either analytically, if pos- 
sible, or numerically. Although the underlying renormal- 
ization scheme is extremely simple the analytical compu- 
tations are sometimes tedious [n |2] . Therefore an alterna- 
tive route to the exact determination of critical exponents 
and other quantities of interest is highly desirable, and 
this is what we are going to present in this letter. In doing 
so we follow a route on which we already traveled suc- 
cessfully for the random transverse Ising chain pl|-fi3|, 
and here we are going to do one step further studying 
random XX and XY models with the help of a straight- 
forward and efficient mapping to random walk problems. 
This mapping is not only a short-cut to the results known 
from analytical RG calculations, it also gives new exact 
results in the off-critical region (the Griffiths-phase |14|) 
and provides a mean to study situations in which the RG 
procedure must fail, as for instance in the case of corre- 
lated disorder [[L5| . Here we confine ourselves to a concise 
presentation of the basic ideas including the determina- 
tion of various exponents for the first time. The technical 
details of the derivations and further results are deferred 
to a subsequent publication Jig ]. 

The model that we consider is a spin-1/2 XY-quantum 
spin chain with L sites and open boundaries, defined by 
the Hamiltonian 



L-\ 



H = J2( J i s iSi+i + J!SfS i 



v \ 
i+i) 



(1) 



where the Sf' y are spin-1/2 operators and the interaction 
strengths or couplings Jf' v > are independent random 
variables modeling quenched disorder. In the case of the 



random XY chain one has two independent distributions 
for the couplings J x and J v , p x and p v , respectively, 
whereas the random dimerized XX-chain has perfectly 
isotropic couplings Jf = Jf = Ji but two independent 
probability distributions for the even and odd couplings 



(i.e. for J 2 i = Jf, and J21-1 = J|i. 

tively. 



_ x ), p e and p°, respec- 



The model @) has a critical point given by [In J x ] av = 
[In J y ] av in the XY case and [In J e ] av = [In J°] av in the 
XX case (here [. . .] av denotes the disorder average). The 
distance from the critical point is conveniently measured 
in the variable 



[lnJ^ e )] av -[lnJ^°)] av 
var[ln J x ( a )] + var[ln J«(°) 



(2) 



where var(x) is the variance of random variable x. At 
the critical point (8 = 0) spatial correlations decay alge- 
braically, for instance in a finite system of length L with 
periodic boundary conditions the bulk-correlations decay 
as 



[C"(Z)], 



[(0|SfS M 



L/2 



|0)]a 



(3) 



for p = x,y,z, (0| denotes the ground state of (Q), 
whereas for a finite system of length L with open bound- 
ary conditions the end-to-end correlations decay with a 
different exponent like 



[C?(L)U = [{0\S?SZ\0)] t 



(4) 



Away from the critical point (8 ^ 0) the infinite sys- 
tem develops long range order. For the XY model 
lim L ^oo[C' i (L)] av = (m») 2 jt 0, with m x > for 8 > 
and m v > for 8 < 0, whereas for the XX model there 
is dimerization for 8 ^ with non-vanishing string or- 
der 1 171]. One can introduce local transverse and lon- 
gitudinal order parameters m XlV and mf also for a finite 
system (with open boundaries) using the off-diagonal ma- 



trix clement 



<•/ Jav 



= [(l|5f |0)] av , where (1| is the lowest 



excited state with a non- vanishing matrix-element p8| . 
Analogous to bulk and end-to-end correlations the bulk 
and surface magnetizations rn^, 2 and m^, respectively, 
behave differently: 



m 



L/2J 



L' 



and 



L" 



(5 = 0) (5) 



where the critical exponents x^ 1 and x^ fulfill the scaling 
relation 2x^ = if and 2x^ = 77^. 

Now we are going to determine the critical exponents 
introduced above. We use the free-fermion representa- 
tion of model (El) to derive exact expressions for the local 
order-parameters whose finite size scaling behavior fol- 
lows then from a mapping to a random walk problem. 

We start with the longitudinal order parameter for the 
random XX chain (of length L with open boundaries; for 
convenience we assume from now on that L is a multiple 
of 4), which is given by |lq| 

L/2-1 k / N 2 

J 2.7-l 



I L/2-1 k 



(-1 k , T v 2' 



fc=lj=l 



21-23-1 



for odd sites. Note that the couplings to the left and 
to the right of the spin that is considered enter this ex- 
pression differently. For the surface (1 = 1), where one 
only has "right" couplings, this expression is similar to 
an analogous result for the random transverse Ising chain 
p] and its scaling properties are related to the survival 
probability of a random walk with L/2 steps. This is 
easy to see for the extreme binary distribution, in which 
J23 = 1 and J2J—1 — A, A -1 with probability 1/2, tak- 
ing the limit A — > (i.e. A -1 — * 00). Due to the occu- 
rance of infinite terms in the sum in the denominator of 
the r.h.s. of ([|) one can easily identify those instances 
that give a non-vanishing surface magnetization: When 
Vfc = 1, . . . , L/2 — 1 : 7 =i J2J-1 < 00 the expression 
on the r.h.s. of (pj) attains a non-vanishing value (typ- 
ically 1 or, less frequently, some fraction 1/n), other- 
wise it is zero. One can represent the disorder config- 
uration Ji, J 3 , J 5 , . . . , Jl-i as one instance of a random 
walk with L/2 — 1 steps by saying that the walker in the 
i-th steps moves downwards if J-xi—i = A and upwards if 
J21-1 = A -1 , as it is sketched in Fig.l. In this way the 
disorder configuration with non-vanishing surface mag- 
netization mf are easily identified: they represent sur- 
viving walks, i.e. walks that never move into the upper 
half. Thus [mf(XX)] av scales like the survival proba- 
bility P surv (L/2) of a random walk with L/2 steps that 
vanishes like L" 1 / 2 , i.e. [mf (XX)] av ~ L" 1 / 2 . Therefore 



x\(XX) 



and ril(XX) = 1 



(7) 



Inspecting the expression (g) for the bulk (I = L/A) 
one sees that now (again for the extreme binary dis- 
tribution) a nonvanishing magnetization TOf,/ 2 _i arises 
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FIG. 1. Sketch of the configuration of odd bonds for a chain 
of length L that gives a non-vanishing longitudinal magneti- 
zation mf ~ 0(1) for the surface spin, i = 1, in (a) and the 
central spin, i — L/2 — 1, in (b). The example is for the 
extreme binary distribution with Jii = 1. Weak couplings 
(J2i-\ — A) correspond to downward steps of the random 
walk, strong couplings (Jii-\ — A -1 ) to upwards steps. The 
walk in (a) has surviving character, it does not enter the up- 
per half plane. In (b) one can identify two random walks each 
starting at the central site, i — L/2 — 1, one to the right and 
one to the left, and each of them has surviving character. 



only if Vfc = 

and Vfc = 1, 



L/4,...,L/2-l : n;=L/ 4 ^-i < 00 

■ • , L/A - 1 : nj=i Jl/2-2j-i < <»■ We 
represent the disorder configuration to the right of the 
central site, Jl/2-i> Jl/2+i> ■ ■ ■ 1 Jl-1, as a random walk 
with L/A steps in the way as for the surface spin. The 
disorder configuration Jl/2-i, Jl/2-3, ■ ■ - > J\ to the left 
is represented as a second (independent) random walk 
also with L/A steps, now counting backwards and with 
the step-direction reversed (i.e. downwards for J = A -1 
and upwards for J = A), since now strong bonds on odd 
sites imply weak coupling of the central spin. For illus- 
tration this representation is depicted in Fig. 1. Now, 
for the bulk magnetization rn z L i 2 _ 1 to be non- vanishing, 
both halfs of the coupling configuration have to repre- 
sent surviving random walks. Thus the probability for a 
non-vanishing magnetization m z L , 2 _ 1 is just the product 
of two survival probabilities (since both walks are inde- 
pendent), i.e. [m| ;/2 _ 1 ]av ~ {-P S urv(£/4)} 2 ~ L^ 1 and 
therefore 



x z (XX) = l and rf(XX) = 2 



(8) 



For the XY chain one has the following exact relation 
|16| for the disorder averaged longitudinal magnetization 

[mf(XT)] av = [{mf(XX)} 1 / 2 ] 2 v , (9) 

which yields immediately the surface and bulk exponents 



for longitudinal order parameter and correlations: Since 
{mf (XY)} 1 / 2 has a non-vanishing value if and only if 
mf(XX) is non-vanishing, one obtains [mf(XY~)] av ~ 
{-Psurv(-^/4)} 2 ~ L" 1 , which means 



c$(XY) = l and r)f (XY) = 2 



and further, 
which means 



b L/2- 



- X {XY% 



{P SUYV (L/A)} 4 



x z (XY) = 2 and rf{XY)=A. 



(10) 
(11) 



For the transverse surface order parameter mf one has 
L6[ the exact formula (valid for the XX and XY case) 




-1/2 



(12) 



which is similar to (g) for mf (XX), apart from the power 
1/2 on the r.h.s., which again is only non- vanishing for 
disorder configurations that represent surviving walks a 
la Fig. la. This then yields without any effort 

x x {XX,XY) = l and rff (XX, XY) = 2 , (13) 

For the transverse bulk order parameter in the XY 
chain we use the fact that the model can be mapped 
onto two transverse Ising models (TIM), with uncorre- 
cted disorder in both chains Pjl6|. Through this map- 
ping one obtains for the transverse correlation function 



c. 



2i,2i+2r 



(0\S5iS5 i+2r \0) the following identity |f 



[C^2 l+ 2r(XY)] av = 4[C* i+r (TJM free )] av 
•[C? i+r (T/M fixed )] av 



(14) 

-2^ {TIM) 



where fixed and free indicated the boundary conditions. 
Since the correlation function exponent is known exactly 
§ to be r] x (TIM) = (3 - \/5)/2 we have: 

x x (XY) = (3-\/5)/2 and f] x (XY) = 3 - VE . (15) 

For the XX chain the two transverse Ising chains have 
perfectly correlated disorder, which implies that the dis- 
order averaged transverse correlations do not factorize 
into two independent averages as in (|l4j). Therefore, 
for the transverse order parameter exponent in the XX 
case we have to use a different route: The first impor- 
tant observation is that the transverse bulk order pa- 



rameter m 



L/2 



1\S^, 2 \0) attains its maximum value 
1/2 if the central spin is decoupled from the rest of the 
system, i.e. when Jl/2~i = Jl/2 = 0- More generally 
we expect that m x L , 2 ~ 0(1) when it is weakly cou- 
pled to the rest of the system. "Weakly coupled" in 
the case of the extreme binary distribution means that 
the bond configuration to the left and to the right of 
the central spin represent both surviving random walks, 
as exemplified in Fig. lb (this is actually equivalent to 



the (exact) condition for the longitudinal order param- 



eter m z L/2 (XX) to be non- vanishing) . 



'L/2 

dence implies 



l L/2 



(XX)] E 



{P s , 



This correspon- 

,(L/4)} 2 ~ L- 1 



from which one obtains 



x x {XX) = 1 and rf{XX) = 2 . 



(16) 



We verified the strong correlation between weak coupling 
and non-vanishing transverse order parameter numeri- 
cally in the following way: We considered a chain with 
L + l sites and the couplings at both sides of the central 
spin were taken randomly in the form of surviving walk 
character, where we used the binary distribution with 
A = 0.1. For such small value of A the surface order- 
parameter averaged over the surviving walk (sw) configu- 
rations [mi] sw was very close to the maximal value of 1/2. 
Then we calculated numerically the order-parameter at 
the central spin and its average value over surviving walk 
configurations [m x L , 2 \ svl as given in Table I. 



L 


2[mf] sw 


2[ TO L/2] SW 


32 
64 
128 
256 


0.994 
0.991 
0.991 
0.991 


0.764 
0.682 
0.647 
0.577 



TABLE I: Surface and bulk transverse order-parameters av- 
eraged over 50000 surviving walk coupling configurations for 
the binary distribution (A = 0.1). 

As seen in the Table the averaged surface order- 
parameter stays constant for large values of L, whereas 
the bulk order-parameter decreases very slowly, actually 
slower than any power. The data can be nicely fitted 
by [to£/ 2 ]sw ~ (InL) -1 / 2 . Thus we conclude that the 
numerical results confirm the exponents given in (|lq), 
however there are strong logarithmic corrections, which 
imply for the average transverse correlations 



[C x (r% 



; ln -1 (r) 



XX - model 



(17) 



These strong logarithmic corrections render the numeri- 
cal calculation of critical exponents very difficult |17|,|l9[ ■ 
In earlier numerical work using smaller finite systems dis- 
order dependent exponents were reported p9j. We be- 
lieve that these numerical results can be interpreted as 
effective, size-dependent exponents and the asymptotic 
critical behavior is indeed described by Eq. ( J1T|) . 

Away from the critical point the correlation length ex- 
ponent v can be determined by the scaling behavior of 
the longitudinal surface magnetization, i.e. (1q) with 1 = 1, 
which can be inferred from the survival properties of a, 



now biased, random walk: [mf' 1 ' (5, L)\ 



P surv (5,L/2). 



A non- vanishing distance (5, see (||), from the critical 
point means that the disorder configurations can be rep- 
resented by random walk that have a drift either towards 
(5 > 0) or away from (6 < 0) an absorbing wall (take for 
instance the extreme binary distribution, in which weak 



bonds A occur with a probability (1 — 8)/2 and strong 
bonds A -1 occur with a probability (1 + S)/2 and com- 
pare with Fig. la). Recalling the asymptotic properties 
of the survival probability of random walks jll[ one gets 
for S > P surv (<5 > 0,1) - exp(-L/£) with £ - 5~ 2 , 
where the characteristic length scale £ of surviving walks 
corresponds to the average correlation length of the XX 
and XY chains: 



[& 



with 



(18) 



For 5 < the drift away from the adsorbing wall 
yields a finite survival probability even in the infinite 
system, P SU rv(<5 < 0,L/2) oc 5, which implies that 
[r<] av ex \8\-tt, with fi' v (XX,XY) = f3f(XX) = 1 
and P{{XY) = 2. Since z^ = ft/v and v = 2 from |[|) 
one reconfirms the results about the surface magnetiza- 
tion exponents in (0), (|o|) and ( |l3| ) 

The typical correlation length, £ t yp can be inferred 
from the scaling behavior of the typical surface magne- 



tization In mi ~ S J -{ha(J| 7 -_ 1 



ln(^2, (e) )} « SL > whi ch 



gives 



K1 



typ 



6 Vtyp with f typ = 1 



(19) 



The critical and off-critical scaling behavior of the low 
energy exciations and dynamical correlations can be de- 
duced from the formula for the gap e |lq| 



L/2- 

e(L)=mfml_ 1 Jl_ 1 JJ 



J. 



v(°) 

2j-l 
LO. rx(e) 
i=l J 2j 



(20) 



which is analogous to a corresponding formula for the 
gap in the transverse Ising chain Q|. At the criti- 
cal point one observes that lne is a sum of L inde- 
pendently distributed random variables with zero mean 
(since 5 = 0), for which the central limit theorem ap- 
plies. Therefore the probability distribution of gaps 
obeys P(lne) ~ L _1 / 2 p(lne/L~ 1 / 2 ) and one uses scaling 
arguments as in [E0[ to deduce the asymptotic (imagi- 
nary) time dependence of the spin-spin autocorrelation 
function G?(t) = [{0\S? (t)S? (0)|0)] av . 



G»(t) ~ (lnr)- 



'/;, 



(21) 



for the surface (a = 1) and bulk (a =bulk), respectively, 
with the critical exponents ry^ as given above. 

Away from the critical point in the Griffiths-phase pj] 
the gap distribution has still an algebraic tail P(e) ~ 
e — l+i/a (S) ^ w j^ n a dynamical exponent z'{8) that varies 
continuously with the distance from the critical point 8 
and is given by the exact (implicitc) formula Jl2[ 



jx(e) 
Jv(o) 



i/z\sy 



(22) 



The dynamical exponent z'{8) parameterizes all Griffiths- 
McCoy singularities occurring in the Griffiths-phase, e.g. 
the spin-spin autocorrelations decay algebraically as 



Gi(t) 



-l/z'(5) 



(23) 



which gives for the susceptibility x^ ~ y— 1+V-* (5) jj_ 
verging for T — > (T = temperature) if z'{8) > 1. 

To summarize we have shown how to obtain a com- 
plete description of the critical and off-critical singular- 
ities of random XX and XY chains with simple random 
walk arguments using exact formulas arising from the free 
fermion description of these quantum spin models. All 
results for the critical exponents are therefore exact. One 
should note that for the transverse bulk order parameter 
exponent for the XY model we referred to a result for the 
transverse Ising model obtained by a RG calculation |lj] 
and for the same exponent of the XX model we showed 
the existence of strong logarithmic corrections. 
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